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This paper presents symmetry classes of the Hartree-Fock (HF) solutions of spin and 
orbital ordered states in a t2g Hubbard model on a two-dimensional square lattice. Using a 
group theoretical bifurcation theory of the Hartree-Fock equation, we obtained many types of 
broken symmetry solutions which bifurcate from the normal state through one step transition 
in cases of commensurate ordering vectors Qq = (0,0), = {n,-K), — ('""lO) ^'Hd 

— (0,7r). Each broken symmetry state is characterized by the presence of local order 
parameters(LOP) at each lattice site: quadrupole moment Q = (Qii Q12, Q23, Qai), orbital 
angular momentum I = {hjhja), spin density s = (s^,s^,s^), spin quadrupole moment 

= (Qi^) <3i2) Q231 Q31) a-nd spin orbital angular momentum Z"^ — (^1,^2 1^3) where A = 
1,2,3. We performed numerical calculations for some parameter sets. Then we have found 
that many types of non-coUinear magnetic orbital ordered states having LOPiQ"^ and Z"^ can 
be the ground state for these parameter sets. 

§1. Introduction 

Perovskite-type i^TiOs and RMOs^^-^^ {R being a rare-earth ion or Y) and quasi 
two-dimensional ruthenate compounds Ca2-a;Sra;Ru04^^"^^ have gained considerable 
interest because of their plentiful phases including various orbital orders. In a pre- 
vious paper^) (hereafter we cite it as I), in order to study these perovskite-type 
compounds, we have considered a triply degenerate t2g Hubbard model on a three- 
dimensional cubic lattice. 

In the paper we have presented a brief review of the general group theoretical 
bifurcation theory of the Hartree-Fock (HF) equation. By listing axial isotropy 
subgroups of the R-reps (irreducible representation over the real number field) of 
the group Gq of symmetry of the system, we have obta ined various types of non 
magnetic orbital ordered states and collinear0 magnetic 1^ orbital ordered states, 



which bifrurcate from the normal state through one step transition in the cases of 
ordering vectors Qp = (0,0,0) and Qq = (vr, 7r,7r). 

In the present paper, to study quasi two-dimensional ruthenate compounds, we 
consider a three-orbital t2g Hubbard model on a two-dimensional square lattice. 



*' All spins are along the z axis. 

**' The phrase "magnetic" means magnetic due to the spin of electrons. 



typeset using VTpJ^.ds (Ver.0.9) 
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We apply the group theoretical bifurcation theory to this model. 

In this paper we consider broken symmetry states with ordering vectors V point : 
Qq = (0,0), M point : Qi = (vr, tt) and X point : Q2 = (tt, 0) and ^3 = (0,7r) which 
can allow states with double Q. Also we consider magnetic states with non-collinear 
spin structure, which were not treated in I. 

This paper is organized as follows. In §2 we give a model Hamiltonian and 
its symmetry group Gq. In §3 wc present a general HF Hamiltonian with ordering 
vectors Qq, Qi, Q2 and Q3, and define the isotropy subgroup of the HF Hamiltonian. 
We present general formulae for the local order parameters (LOP) at the lattice site 
m: the charge density, the spin density, the quadrupole moment, the orbital angular 
momentum, the spin quadrupole moment, and the spin orbital angular momentum. 
In §4 wc present R-reps of the symmetry group Gq of the system over the HF 
Hamiltonian space. 

In §5 we present symmetry classes of non magnetic orbital ordered states by 
listing axial isotropy subgroups of R-rcps which do not break the spin rotation sym- 
metry. In three examples, we show how to derive the canonical form of the HF 
Hamiltonian and occupied orbitals from the isotropy subgroup of a state. We show 
that all states, which break the time reversal symmetry, have orbital angular mo- 
mentum and almost of all states, which do not break the time reversal symmetry, 
have the quadrupole moment. 

In §6 we present symmetry classes of magnetic orbital ordered states by listing 
axial isotropy sTibgroups of R-reps which break the spin rotation symmetry. From R- 
reps which break time reversal symmetry, we obtain states which have spin densities 
or spin quadrupole moments as LOP. From R-reps which do not break time reversal 
symmetry, we obtain states which have spin orbital angular moment as LOP. 

In two collincar examples, in which hold the spin rotation symmetry around 
the z axis, and the time reversal symmetry is broken, we show how to derive the 
canonical form of occupied orbital and explicit form of the spin density or the spin 
quadrupole moment as LOP. In three non-collinear examples, in which both of the 
spin rotation symmetry around the z axis and the time reversal symmetry are broken, 
we show how to derive the canonical form of occupied orbital and explicit form of 
the non-collinear spin density or the spin quadrupole moment as LOP. 

In §7 we report numerical results for some parameter sets. There we show that 
various non-collinear magnetic states can be the most stable in all states described in 
§5 and §6 for these parameter sets. In §8 a summary and discussion are given. In this 
section we discuss two examples of states bifurcating through two step transitions 
from the normal state. It is shown that these are the coexistent states of the spin 
density and the quadrupole moment. The notation used in this paper follows the 
one in I. 
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§2. The model Hamiltonian and its symmetry 

We consider the three-orbital t2g Hubbard Hamiltonian on the two-dimensional 
square lattice (lattice constant= 1): 

3 3 

= X] X/ X] ^i^nis^nis ~ X] ^fj^liis^{n+a)js ~ ''^nis 

n i=l s n,a i,j=l s n,i,s 

+ U ^ n^iinnil + — ^ ^ ^ f^nis'f^njs' 
n,i n ijLj s,s' 

+ ^1212Y1 "Ls"ijs'ams'arijs + ^J2Y.Y1 «Ls«L'«r^is'ar^i., (2-1) 
n ijLj s,s' n i^j s,s' 

where a^^^ is the creation operator for a t2g electron with spin s in the i-th (i = 
1 for dyz,i = 2 for d^x and i = 3 for dxy) orbital at site n and n^^^ = aJij^a^js, Ci is 
the crystal field for the i-th orbital, a is the vector connecting nearest neighbor sites 
, that is, a = (1, 0), (0, 1), (— 1, 0), (0, — 1) and t° is the nearest neighbor hopping 
integral between i and j orbitals along a direction, U is the intra-orbital Coulomb 
interaction, U' is the inter-orbital Coulomb interaction, J is the exchange integral 
and J' is the pair hopping interaction. U = U' + J + J' is derived from rotational 
invariance in orbital space and J = J' from the evaluation of Coulomb integrals. 
Three t2g atomic orbitals (pi{r), 4'2{f) and (j)z{r) of an atom located at the origin 
are defined by 

(t)i{r) = dyz = f{r)yz, 
(p2{r) = dzx = f{r)zx, 

Mr) = d,y = f{r)xy, (2-2) 

where r = {x,y,z), r = \r\ and /(r) is a spherical symmetric function. Three t2g 
atomic orbitals of an atom located at a site n are given by (j)i{r — n), 4'2{f — n) and 
03(r - n). 

The symmetry group Gq of the Hamiltonian 7i in (j2-l|) is given by 

Go = P X 5 X R, (2-3) 

where P = L(ei,e2) A D^h is the space group of the square lattice, I/(ei,e2) is a 
two-dimensional translation group with basis vectors ei and e2, and A denotes the 
semidirect product. D/^h is the point group of the square lattice. S is the group of 
spin rotation, R = {E, i} is the group of time reversal. 

Three orbitals (piir) , 4'2{f) and (t)j,{r) have following symmetry properties for 
P G Dih- 

2 

p . Mr) = Mp'' ■r) = Yl (t>j{r)D^f'\p), for i = 1, 2 

p ■ Mr) = Mp~^ ■ r) = X^''''\p)Mr), (24) 
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Table I. R-rep matrices D^''^ of Dih 
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V -1 7 
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-1 






C2z 
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C2x 
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-1 
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-1 




C2y 


1 


-1 
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-1 




C2a 
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-1 


-1 
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(° o) 


C2b 


1 


-1 




1 


o) 




(p), for? 


5 G £>4, 


I : inversion. 





where D\f'\p){x^^^'Hp)) 

are R-rep matrices of Eg{B2g). The equations 
be expressed in a more compact form 

3 

P- (t}i{r) = ^(l}j{r)Dji{p), i = 1,2,2, 



where 



D{p) 



d[^^\p) D'i^\p) 



\ x^^^'\p) 

Action of p G Dah on (pi [r — n) are given by 

p-(l)i{r -n) = (pi{p'^ .r-n) = (l}i{p~^ ■ {r - p ■ n)) 

3 

= ^'Pi{'>' -p-'<^)Dji{p)- 
Then actions of p € -D4/1 on {a^nis^'^nis) given by 
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Actions of translation T{n) £ L{ei,e2) with a vector n = nie\ + 71262(711 and 
n2 are integers) on {a^rnis-:'^mis} ^"^^ given by 

-"L. = «(r.+m)i.' • = (2-9) 

Actions of spin rotation u{n, 9) G S, time reversal t G R on {aj^js, Ornis} defined 

byl2),14) 

u{n,e) ■ a^^^^ = ^u(n,e)s'sa^rnis'^ u{n,e) ■ a^^, = ^{«(n, ^)y,}*a^j,,, 

s' s' 
^ • = ^*«Lt' * ■ = 2*«miT' (2-10) 

where 7i(n, 0) is a spin rotation by around the n axis and is given by 2 x 2 unitary 
matrix; 

6 6 
u{n, e) = cos(-)l2 - i{cr ■ n) sin(-), 

<T = ((7^, cr^, (T^) : Pauli matrices, 

I2 = 2 X 2 unit matrix, (2-11) 

and z is a complex number and z* is a complex conjugate of z. 

Prom invariance of the Hamiltonian Ti,, we obtain following conditions for 
ei and i"- 



ei 


= £2 = (5, 


€3 = 








^(1-0) 

''ii 


_ /(-i,o) 


_ .(0,1) _ 
— '-22 — 


.(0,- 

'-22 


■1) 


= ii, 


^(1,0) 
'-22 


_ /(-i,o) 

— '-22 


_ .(0,1) _ 

— I;^^ — 


.(0,- 

''11 


■1) 


= i2, 


.(1,0) 

"-as 


_ .(-1,0) 
~ "-as 


_ .(0,1) _ 
~ ''33 ~ 


''33 


-1) 


= h, 




= 0, in 


other cases 







(2-12) 



where 5 is the level splitting between 03 and (4>i,4>2) orbitals. 
Using Pourier transformations 



'4ks - XI ^''^ ""«is' "'ifes - X ^ ''^ ""ms' (2-13) 



we obtain 7^ in the momentum representation 

3 



k.s i=l 



1 ^ 

fc,fe',q i,j,n,m=l s,s' 



6 



M. Hamada, A. Nakanisi, A. Goto and M. Ozaki 



where is the number of lattice sites and k runs in the first Brillouin zone: 
{k = {ki, ^2) I — TT < ki, k2 < vr}, and 

tii{k) and {i{k + q)s,nk's' \ V \ jks^ra{k' + q)s') are given by 

tii(fc) = —5 + 2(ti cos ki + t2 cos k2)., 
t22{k) = —5 + 2{t2 cos ki + ti cos ^2), 

t3s{k) = 2t3(cos ki + cos ^2), (2-15) 



{i{k + q)s, ik's' \ V \ iks, i{k' + q)s) 



U 

iV' 
U' 



{i{k + q)s,jk's' I V I iks,j{k' + = — , i ^ j, 

{i{k + q)s,jk's' I 1/ I ifcs, + q)s') = i / j, 

J' 

(i(fc + q)s,ifcV I y I jks,j{k' + g)s') = — , j, 

{i{k + q)s,nk' s' \ V \ jks,m{k' + q)s') = 0, otherwise. (2'16) 

From (gSD, (gl]), (pTUjl and OTHD . actions of p G D4h,T{n) G ^(61,62), 
u{n,6) e S,t eK on {alfc^'^ifes} ^'^^ given by 

3 3 

Tin) ■ = e"""-"aL' ^l'^) " ^iks = e^'^'Xfcs, 

2 2 

u(n, 6*) • = ^ u(n, e)s's4ks' ' ^) ' ^iks = X] ^)s'^"ifc^' ' 

s'=l s'=l 

§3. Hartree-Fock Hamiltonian and its isotropy subgroup 

In this paper we consider HF solutions with four types of ordering vectors Qq = 
(0,0), Qi = (vr,7r),Q2 = ('^)0),Q3 = (0, vr). Thus a general HF Hamiltonian Hm is 
written as 

3 3 3 2 

Hm = HK + J2J2J2Y.Yl 4(^)«WQ,)s^ss'«,fcs', (3-1) 
1=0 i,j=l X=0 s,s'=l k 

where Hk is the kinetic energy written as 

3 2 

^x = ^^^{-in(fc)-/i}aL«ifcs- (3-2) 

fc 2=1 S=l 
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From the Hermite condition of Hm we have 

x'^ik + QiY =xf,{k). (3-3) 

,, . . f fc^ in (\^- 1 h sntisfiA fnllnwincr 



x'^(fc) in (|3-ip satisfy following SCF conditions 



X, 



(^) = E E WL^^{k,k')pt^{k'), 1 = 0,1,2,3 

k' m,n=l 
3 

(^) = E E Yl,^^{k,k')p'^Jk'), 1 = 0,1,2,3, A = 1,2, 3 

fc' m,n=l 



(3-4) 



where 



iA 1 t A 

s,s'=l 



,„j^(fc, fc') = 2(i(fc + Qi)s, nk's I y I jks, m{k' + Q,)s) 
— {i{k + Q/)s, nfc's | V \ m{k' + Qi)s,jks) 
Yl^j„,{k,k') = -{i{k + Qi)s,nk's I 1/ I m{k' + Qi)s,jks) 



(3-5) 



where (yl) denotes the expectation value of A in the ground state of the HF Hamil- 
tonian Hm and p^^{k) satisfies 



p%ky = p[>;{k + Q 



3 



l)^ 



(3-6) 



A=0 



From (I2-16|) we have for i ^ j 

U 



WUk,k') 


= Wan = 


Wlji,{k,k') 


= ^^ijij ~ 


Wl,,{k,k') 




Wi,^,{k,k') 


= Wijji -- 



iV' 

2U' - J 

N ' 

f 

2J-U' 

N ' 



^iiiiik, k ) = Yiiii - 

Y^j^j{k, k ) = Yijij 
Y^j^jAk, k ) = Yiijj 



^ijjiik, k) — Yi 



U 

'iV' 

J_ 

"iV' 

,r_ 



YinjTni^i = Yinjni = 0, otherwise (3-7) 



From (|3-4|) and (|3-7|) we see that xMk) are independent of k and given by 



•^ij{k) — Xj^j — A'^ ^ ^ ^injmRmm ^ — Oj 2, 3 



m,n=l 
3 



x'^ik) = xfj = N YinjmR'L, 1 = 0,1,2,3, A = 1,2, 3 (3-8) 

m,n=l 
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where 

4' = 4E'^''(^)- (3-9) 



N 

k 



We denote a 3 x 3 matrix whose component is x'^ [Rlj) by x^^{R^^). From 



(IMD . (13^ . (13^ . we see that x^^ and i?'-^ are Hermite matrices. 
From (|3-8|) the HF Hamiltonian in (j3-l|) can be written as 

3 3 3 2 

^- = ^^ + EEEEE ^'M(fc+Q,).^.' (3-10) 

i=0 jj=l A=Os,s'=l fc 

From (13-lOp and the Hermiticity of x''^, Hm is characterized by 3 x 3 Hermitian 
matrices a;'^ (16 matrices in all, corresponding / = 0, 1, 2, 3, A = 0, 1, 2, 3). Thus Hm 
is specified by a vector in the HF Hamiltonian space Whf over real number field 
R: 

WnF = I E(4(fc+Q,)s"jfcs' + a]fcs'ai(fc+Q,) J' ^i^4ik+Q,)s''jks' " i^ks'^^ik+Q^^^ \ ■ 

K k k J R 

(3-11) 

where / = 0, 1, 2, 3, i, j = 1,2,3 and {^4, B, - ■ ■ , }r, denotes a vector space with bases 
A,B, - ■ ■ over the real number field. 

The HF energy is expressed in terms of i?''*' 

3 

^FH = (^^k)' + Y1 W.njmRfiRian 

i,n,j,m=l 
3 3 

i,n,j,m=l A=l 
3 

{HKy = -2Y,Y.^n{k)p'^iik). (3-12) 

fc i=l 

The SCF condition (j3-8p corresponds to the extremum condition of Epii}^^ 
Actions 5 G Go on Hm are defined by 

9 ■ = Hk 

3 3 3 2 

+EEE E E(4')^*H^^-Wq,)J(^-')^*^(^ ^3-13) 

1=0 i,j=l A=0 s,s'=l k 

where A^*^ denotes the complex conjugation of a complex number A in the case of 
anti-unitary g which contains time reversal t. Note that 

3 2 

5 • = E E E i-^-^'^) - • 4ks)i9 ■ = Hk. (3-14) 

k i=l s=l 
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We define the isotropy subgroup G{Hm) of Hm by 

G{Hm) = {geGo\9-Hm = Hm}. (3-15) 
Actions of 5 G Go on Rf^ are defined by 

i9 ■ R%j = ^ E E ((9-' ■ a]^k+Q0s)(9-' ■ aife.'))(*V,\,. (3-16) 

k s,s'=l 

In previous papers, -^^^ we have shown that for g G G{Hm) 

{(9 ■ aljig ■ a.^,^,))(*^ = {al,a^^.^,). (3-17) 
Then we obtain for g G G{Hm) 

g-R^^ = R^^. (3-18) 

For subsequent uses we hst exphcit forms of Go actions. For p G i)4h, T{m) G 
LQ,u{n,9) G 5,t G R 



p- 


R'^ 


= D{p)R^ ''^^^D^{p), 


T{m)- 


R'^ 




u{n, 9) 


■R'° 




u{n, 9) ■ 


R'^ 


3 

= Y.R{n,-9)yxR'^'^ 






A'=l 


t 






t- 




= -{R'^)*, A = 1,2, 3, 



(3-19) 

where {p~^ ■ I) is defined such that 

Qp-i.l = p-' -Ql, (3-20) 

and R{n,9) is the rotation matrix by 9 radian around n = (ni, 712,713) axis in the 
three dimensional EucUd space and is given 

byl9),20) 

R{u{n,9)) =R{n,-9) = 

cos^ + (1 — cos0)nf (1 — cos^)niri2 — n3sin^ (1 — cos0)7Zi7Z3 + n2sin0 
(1 — cos ^)ni7i2 + J^s sin0 cos + (1 — cos 0)n2 (1 — cos0)n2n3 — ni sin0 
(1 — cos^)nin3 — ri2sin^ (1 — cos^)ri2ri3 + nisin^ cos ^ + (1 — cos 6')?23 

(3-21) 

We define density matrices D**'(m) = {Df?' {i,j = 1, 2, 3, s, s' = 1, 2)} at a site 
m as follows: 



Dtfim) = {a^^.^,a mis)- (3-22) 
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Since for all states with ordering vectors Qi {I = 0, 1, 2, 3) 

(alk's'^iks) = 0, for k'^k + Qi, (3-23) 

using ()3-6p we obtain 

fc,fc' 



k 1=0 

k 1=0 A=0 

3 3 

^^e-Qr-^A,^a_ (3.24) 



1=0 \=o 



Thus we obtain 



3 3 



D''\m) = ^^e-'^^-"'ai,R^\ (3-25) 

1=0 X=0 

Then we obtain exphcit expression of density matrices as follows: 

3 

1=0 
3 

1=0 
3 

1=0 
3 

D^^m) = ^ e"^^'-"^(i?'^ + iR^^). (3-26) 



i=0 



From (I3T8P we can see that symmetry properties of density matrices D^^ (m) are 
determined by G[Hm)- 
Using notations 

A\m) = (A\{m),A\{m),A\{rn),A\{rn),Al{rn),Al{rri) 



— ('^Llf '^Tri2t' "rn,3T' ^TTilJ,' '^m2J.' '^Tri3i) ' (3-27) 

we define generalized density matrix T)(m) whose component is given by 

V^j{m) = U]{m)Ai{m)) . (3-28) 
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From (I3-22I) we obtain 

1^{rn)-^ ^^^^^^ ^^^^^ j. (3 29) 

The 6x6 Hermitian matrix X>(m) is diagonalized by a unitary matrix U as follows: 




















\ 





As 






















As 






















A4 






















As 







V 














Ae 


/ 



Defining ^^(m) = (s|(m), 5|(m), 5|(m), ^^(m), 5|(m), 5^(m)) by 
B^{m) = A\m)U, B{m) = U'^A{m), 

we obtain 



(3-30) 



(3-31^ 



{B\{m)Bk{m)) = ^ {A\{m)UuUlA,{m)) 

6 

= Y,mA{m)A,{m))ul 

6 

= ^ ulp{m)jiUa = 6kiXi. 



(3-32) 



Thus occupied spin orbitals -0; and their occupation numbers n^ipi) are given by 

3 



i=l 



n{il)i) = A;. 



(3-33) 



where Z = 1, . . . , 6. 

Here we present formulae for the local order parameter(LOP) at a site m. The 
charge density at a site m: d{m) is expressed by 



j = l ss' 



(3-34) 



The A th component of the spin density at the m: s (m) is written as 

3 



\m) = IYj J2^aL..a^,.>)(Ti,,X = 1, 2, 3. 



2 / > / j \-^mis""mis' 1"^ ss' 
j=l ss' 



(3-35) 
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The component of the quadrupole moment at the site m: Qij{m) are written 
as 

2 

Qii(m) = /i (2(aLisarriis') - (aL2sarri2s') - (^LssamSs')) o-L' 

s,s'=l 
2 

Q22{m)=h Y (-(aLlsarril^') + 2(a|;„2sarri2s') - (aLss^rriSs')) ^L' 
s,s'=l 
2 

<333(m) = (-(oLl^OrriU') " {(^'^rn.2s'^m.2s') + 2(0^.3^ a^g^, )) (T°^, 

s,s'=l 

Qij{m) = h X](«L.«rriis' + aLi.arri».')^L'' « 7^ J> (3-36) 

ss' 

where 

„2 „ 2 



h = j dr4>i{r){x -y )<pi{r), 

l2 = 3 j drMr)xyMr). (3-37) 

The derivation of (I3-36P is given in Appendix A. 

The component of the spin- quadrupole moment at the site m: Q^-{m) are 
defined by, for A = 1, 2, 3, 

2 



TTl3s'/ ; '^SSM 



Qll(^) = \h (2(aLl.arril.') - (aL2s«rri2s') " (^Ls^^^ 

s,s'=l 

^ 2 

Q22("^) = Y (-(^LlsOrrils') + '^{a^rr^2s(^m2s') " (aL3sarri3s')) 

s,s'=l 
^ 2 



2 

s,s'=l 

Qi'il"^) = y X]^"^riis«mjs' + aLjs«mis')f^ss'> ^ ^ J, A = 1, 2, 3. (3-38) 

ss' 

This type of order parameter has been treated in a paper by Shiina, Nishitani and 
Shiba^^^ as the coupled orbital and spin morment in the case of the superexchange 
model of Cg orbital. 

In our system with tetragonal symmetry we use Q|(m) and Q2'^{m) instead of 

Qii and Q^j(i = 1, 2, 3), where 

^ 2 

Qli-m) = -(Qll - Q22) = ((«tils«rrils'> - (aL2sarri2s'>) ^^ss'' 

ss'=l 

1 I ^ 

Ql\rn) = -(Q^i - Q^2) = y E ((«Ls«rr^ls') - («L2s«rr^2s')) ^ss'- (3-39) 

ss'=l 
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The i-th component(i = 1, 2, 3) of the orbital angular momentum at the site m: 
li{m) are written as 

h{m) = {i{a}m.2s(^rruis') " ^'\o)rru^s'^'m2s')) 

ss' 

him) = ^ {i{a^m3s0'7nls') - ^(«Llsam3.')) 
ss' 

him) = {i{al,uam2s') - ^(«L2samis')) ^L'- (3-40) 

ss' 

The derivation of (j3-40p is given in Appendix B. 

The i-th component(i = 1,2,3) of the spin-orbital angular momentum at the 
site m: li{m) are defined by, for A = 1, 2, 3, 

^l{m) = (^(«L2sOrri3s') " Ko^mZs(^-m2s')) ^^ssM 

ss' 

^2 ("^) = \Y^ («(«L3.arrils') " iio^mls^^-mZs')) ^^^s'^ 
ss' 

^3 = \Y^ {Ko^mls<^-m.2s') " Ka^m.2s(^-mls')) ^^^s' ■ (3-41) 

ss' 

Here we express these local order parameters in terms of Fl}^. Using (|3-22p and 
()3-25p we obtain 

3 

E^«Us«-s')^s.' = 2 e- '''^''^Rfr (3-42) 

ss' 1=0 

From (|3-42p . we obtain 

d{m) = 2YYe-^^rrr.j^j] 

1=0 j=l 

s\m) = Yj:e-^^^-R% 

1=0 j=l 

Ql{m) = 2hYe-^^^-^ (<i - 4^2) 
1=0 
3 

Q,,(m) = 2/2 Y (4° + ^f) ' ^ ^ 

1=0 
3 

Qf{m) = hY (^n - ^22) 

1=0 
3 

1=0 
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1=0 1=0 
1=0 1=0 

k{m) = 2^e-Qr-(-z)«2 - Ikm) = ^e-«^-(-i)« " 

«=0 1=0 

(3-43) 

Here we give the physical meanings of Qfj and Ij A = 1, 2, 3). As an example 

we consider the case of and if. The quadrupole moment Q\2{m){Q[2{'fn)) by 
the up spin (down spin) electrons at the site m is written as 



Thus we obtain 



Q\2^m) = h (aLii«rri2i + aL2i«rriii ) • (3-44) 



Qi2(m) = Q\2{m) + Qial"^), 

Q?2(rn) = \ (QUm) - Qi2{m)) . (3-45) 



Then we obtain 



QU^) = \ (Qi2(m) + 2Q\2{m)) , 

Qi2i^) = \ {Quim) - 2Qi2{m)) . (3-46) 

From (j3-46p the existence of (5x2 ("^) implies the different quadrupole moment for 
different spin. 

The orbital angular momentum l\{m){l[{m)) by the up spin (down spin) elec- 
trons at the site m is written as 

/j(m) = i <^a|^2T"rri3T ~ "L3T"rri2T 

liim) = i (^al^2i<^m3i - al^3ia^2i) ■ (3-47) 
Thus we obtain 

ll{m) = ^{h{m) + 2ll{m)), 

li{m) = ^{h{m)-2ll{m)). (3-48) 

Thus the existence of lf{m) implies the different orbital angular momentum for 
different spin. We note that other Q^j{m) and lj{m) have physical meanings similar 
to the above cases. 
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§4. R-reps of Go in Wuf 

First wc give some notations and definitions in reference to tlic group theory. 
We denote an R-rep of a group G as where 7 labels an R-rep. Let (Pig) be the 
R-rep matrix of G^ corresponding to g £ G and W{G'^) be the representation space 
of G'' spanned by {ij,^^ - ■ ■ , /«) over the real number field: 

w{G^) = {q,il... JZU. (4-1) 

Then for 5 e G 

n 

g.q = Y,d^{g)al. (4-2) 

i'=l 

Using real numbers Xi{i = 1, 2, • • • , n), a vector v G W(G'^) is written as follows 

v = Y.xdl (4-3) 

i=l 

The isotropy subgroup ^{v) of v e W{G'') is 

Eiv) = {geG\g-v = v}. (4-4) 

Let i7 C Go be a subgroup of Gq. The fixed-point subspace of H in W{G'^) is 

Fix^(ii") = {ve VF(G^) \g-v = vforallge H}. (4-5) 

An isotropy subgroup with one-dimensional fixed-point subspace is called an axial 
isotropy subgroup }'^^ 

According to group theoretical analysis of the HF equation, ""^^^ instabil- 
ities of a HF solution (characterized by HF Hamiltonian H^) is labeled by R-reps 
of the isotropy subgroup of Hm- 

Using the equivariant branching lemma^^^ in the group theoretical bifurcation 
theory/^)' we can show^^' that if an instability of a HF solution characterized 

by an R-rep Gq occurs, there always exists a branch of a HF solution which bifurcates 
through the instability and has the axial isotropic subgroup in W{Gq). Thus we can 
enumerate broken symmetry states bifurcating from the normal paramagnetic state 
by listing axial isotropy subgroups of each R-rep Gq of Gq. 

Since we consider HF solutions with three types of ordring vectors, F point: 
Qq = (0,0), M point: Qi{it,tt) and X point: {Q2 = {t^,0),Q^ = (0, vr)}, we present 
R-reps of Gq in the representation space Wuf with ordering vectors Q^, (i = 0,1,2,3). 
An R-rep of Gq = P x S x R in Whf is written as Kronecker products of R-reps of 
P, S and R as follows. 



Go = -P <8) 'S' (g) R. 



(4-6) 
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Thus relevant R-reps of Go in Whf are written as 

(4-7) 

where is the identity representation and is a three dimensional R-rep of S 
written as 

S^^\u{n,e)) = R{u{n,e)), (4-8) 

where R{u{n,6)) is the three-dimensional rotation matrix defined in p-2ip . R*^'''' is 
the identity representation and R^^^ is a one-dimensional representation such that 

R(i)(S) = 1, RW(t) = -1. (4-9) 

p(^i)^p(A^i)-j g^^g R-reps of P with ordering vector Qq = (0, 0)(Qi = (tt, vr)) and j is 
the label of R-rep of the little co-group D^h of Qq{Qi), that is , j = Aig, A2g, Big, A2g, 

Eg0. The Rep p^^'^^ are R-reps of P with ordering vector Q2 = (tt, 0) and 7 is the 
label of the R-rep of the little co-group D2h of Q2, that is, 7 = Ag, Big, B2g, Bsg. 
The R-rep matrices of P^^^\ p^^^^^ and p^^^^ are given in Table ITTl 

In Table Im] and Table HV] we hst bases in Wuf of Gf^'''''"\ gJ^^^'^'^^ and 

gI)^^'^'''\ Then bases h{A, ^i,u)m,x,i of P^"^^ S^''^ R(^) transform as follows 
for p e Dih,T{m) G Lq, u{n, 0) G S and r G R 

pT{m)u{n, e)r ■ h{A, fi, v)m,x,i 

= E E^-'-(p^(^))^A'l(^(^>^))R^'^HoMA/^,^)m',A',i, (4-10) 

m'=l A'=l 

where | A \ denotes the dimension of a R-rep A. 

In the following sections we show that R-reps G^''^'^\x 7^ rAig,MAig,XAg) 
derive states with quadrupole moment, R-reps G^q^'^'^\ in which the time reversal 
symmetry is broken, derive states with orbital angular momentum, R-reps 
A 7^ rAig,MAig,XAg) derlve states with spin quadrupole moment, and R-reps 

g'^'^'^\ which hold the time reversal symmetry, derive states with spin orbital an- 
gular momentum. 

§5. Symmetry classes of non magnetic orbital ordered states. 

In this section we consider broken symmetry states derived from R-reps g\^^'^''^\ 
Ci''^'''^'^ and which hold spin rotation symmetry S. In order to list non 



*•* The HF Hamiltonian space Whf is spanned by only gerade(even) bases. 
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Table II. R-rep matrices of P, S and R 



group R — rep matrix 



P P^'^'\pT{m)) = D^J\p) peD^h 

p'''""\pT(m)) = D^J\p)e-'Q^-'^ p G D^h 

" \ X^"'\C2apC2a)e-^^^-'^ J P 2h 



S S'™(^t(n,0)) 


1 


u(n,6i) e S 


s'-''\u{n,e)) 


R{u{n,9)) 




R R(°'(t) 


1 


t £ R 


R(i)(*) 


-1 




(l)T(m)eio = i(ei,e2) 






(2)_D'--'' = an R-rep matrix of Z)4h 






(3)x''''' (p) = an R-rep matrix of D 


2h 




{4:)R(u{n,6)) = the matrix defined 


in (|3-21|) 





magnetic ordered states with ordering vectors Qi (/ = 0, 1, 2, 3) bifurcating from the 
normal state, we present the axial isotropy subgroups of the R-reps q'^^'^''^^ ^Q^^^^'^''^^ 
and G'^"''^''^\ The axial isotropy subgroups of R-reps G^f'''^''^^ and Gq^^'^''^''^^ are 

listed in Table El The axial isotropy subgroups are listed in Table IVTl 

All isotropy subgroups in Table IVl and IVII contain S. Thus from S invariance of 
R^>', only are non-zero. From (j3-26p the density matrix is given by 



D{m) = D^\m) = D^^{m) ="^e-'^'-'^R^°. (5-1) 

;=o 

Since the R^^ is an Hermitian matrix, D{m) is diagonalized by a unitary matrix U 

as 

U^D{m)U = A, 

/ Ai \ 
A = A2 . (5-2) 
V A3 / 

We define 

-^\ns — ('^Irils' '^Tn2si '^Tn,3s)' 

\ns — ("^Inlsi '^Tn2si '^Tn2s) = ^Ins^- (^'3) 



a 
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Table III. Bases of G^^^'"'''' and G^,^^'"'") in VKhf 



R — rep Bases in Whf 



f = EfcE.s'(«2fcaSfc.' +«L.Sfc.')^. 



A 



° 1 ft(rSg,/i,M)2,„(A),l =EfeE,.'Kfe,a3fc.' +«3fe.'*lfc.')°'«^ 



;5(ri5„M,l-M) / h(rEg,H, 1 - M)l,n{A),l = EfcE.s' *(«2fc.Sfc.' - "3fc3"2fc«')^»s' 

h{rEg,ij., 1 - M)2,n(A),i = EfeEa' ^(alfe^Sfey - 4fes«ifey)<^^^' 

/l(M^?g,M,M)l,n(A),l = EfeE..' «3(fe+Qi.)Sfe.''^»V 
^(MA,„M.1-M) /l(M^2g,M,l-M)l,n(A),l =EfcEw«(aI(fc+Q^,)"2fcs' -«2(fc+Q,s)"lfc.')f^s 

^(MSi,.;x.;x) = EfcE..'(«I(fc+Q^,)aife,, - 4(fe+Qi,)«2fcs')'^.^' 

^(MB,„^,^) hl^MB2g,fl,n)l,n(X),l = EfeE.»'(«I(fe+Q^,)a2fe.' +«2(fe+Q,.)«lfe.')<^-' 



^(M£„^,p) f h{MEg,H, M)l,n(A),l = Efc E..' (4(fc+Q,»)"3fc.' + «3{fc+Qi»)Sfe.')f^.\' 

° '»(^^S.M,M)2,n(A),l =EfeE.y(«l(fe+Q^,)a3fe.'+4(fe+Qi»)«lfe.')'^-' 



^(Mi5„;..l-M) f 1 - M)l,n(A),l = EfcE.s' «(«2(fc+Qi .)Sfe.' " «3(fc+Qi .)«2fe.' 

A = for 11 = 0, and A = 1,2,3 for ^ = 1. 

n(A) = 1 for A = 0, and n(A) = A for A = 1, 2, 3. 



These mean 



i'=i 

3 

"mis = X ^rnj'sUpj- (5-4) 



j'=l 



Thus we obtain 
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Table IV. Bases of G, 



(X7,M,i') 




in Whf 



R — rep 



bases in W^hf 



h{XAl,fX, fxh,nWA = Efc Ess' al(k+Q^)s''2ks'<'^ss' 
h{XAl,,l, n)l,n(X),l = Efc E..' 4(fc+Q,),«2fc,''^«^' 

_ hiXAl, m)2.„(a),i = Efc E..' «I(fe+g3),«ife.'^.^' 

h{XAl,H,n)i^ri^X),l = Efc Ess' 4(fc+Q2)s'*3fcs'^ss' 

ll(XAl.fJ. fj)2.„(\).l = Efc E,.,,' "!„fc I Q„;J':;k.'"'^-^' 



h{XBlg,IJ.,IJ,)l^^(^XU = Efc Ess' K(fc+Q,)s«2fes' +«i(fc+Q,)s«lfcs')«^=^ 

^(X Big, p., p.) 

/l(XBig,/X,/x)2,„(A),l =EfcEss' («2(fc+Q3)s'^lfe.' +"I(fc+Q3)s"2fcs')^s^' 
h{XBlg, H, 1 - m)i.„(A),1 = Efc Ess' »K(fe+Q,)sSfcs' - "2(fc+Q,)s«lfe»')'^- 

/s(X-Bi„,(i,l-;i) ' 

h{XBl„H, 1 - M)2,n(A),l = Efc Ess' «(«2(fc+Q3)s«lfc.' " '^l(k+Q^)s'^2ks') ''s 
/l(XS29,M,M)l,n(A),l =EfcEss' («I(fc+Q2)sSfcs' +4(fe+Q2)s'*lfes')f^»^' 

^(XB2g,M,M) f 

_ /l(XS2s, /X, M)2,n(A).l = - Efc Ess' («2(fc+Q3)sSfcs' + "■l(k+Q,)s''2ks')'^s 



^(XB2g,M.l-M) 



^(XB3g,(i,l-(i) 



h{XB2g,H, 1 - M)l,n(A),l = Efc Ess' «(«I(fe+Q,)s«3fcs' " ^(fe+Q^js^lfe J 
HXB^^H, 1 - M)2,n(A),l = - Efc Ess' «(«2{fc+Q3)»"3fes' " «3(fc+Q3)^"2fe-')''- 
h{XB3g,H,n)l,n(X),l = Efc Ess' («2(fc+Q2)sSfcs +4(fc+Q2)s"2fcs')'^^^' 
h{XB3g,H, M)2,n(A),l = - Efc Ess' («I(fe+Q3)s"3fes' +«3(fe+Q3s'*lfe.')'^^''' 
h{XB3g,H, 1 - M)l,n(A),l = Efc Ess' *(«2(fc+Q2)sSfcs' " "U+Qs)-"^*'-')''^^''' 
/l(XS3s,M,l-M)2.n(A),l = -EfcEss'«(«I(fe+Q3)s"3fe.' -4(fc+Q3s'*lfe.')^- 



A = for /i = 0, and A = 1, 2, 3 for /x = 1 

n(A) = 1 for A = 0, and n(A) = A for A = 1, 2, 3. 
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Table V. Axial isotropy subgroup and its Fixed point subspacc of G,,^^'*'''^^ and Gq 



R-rep 


axial isotropy subgroup 


Fixed point subspace 


^(rAig,o,o) 


G(r^ig,0,0) = 


= ID4fi,LoST{, 


Sh \ 
i,'ll,l,lJ-R 


^(/'^2s.0,l) 


G{rA2g,G,i) -- 


= MxC ihLoS 


4/7,1 1 1 Vt3 


^0 


G(rBig,0,0) = 


= D2hLoSTl, 


f hi 1 1 1-R 


^0 


G(rB2g,Q,Q) -- 


~ D2ahLoSI{. 


l/ll 1 1 Ir 


^0 


G(r£g, 0,0)1 : 


= C2xhLoSIl 


{^1.1, i}r 




G(rs9, 0,0)2 = 


= C2ahLoSI{. 


{^1,1,1 ~ ^''2,1,i}r 




G{rEg,o, 1)1 : 


= MzC2xhLoS 


i h 1 1 1 




G(r£<„ 0,1)2 = 


= MzC2ahLoS'R, 


{^1,1,1 ~ ^^'2,1,i}r 




G(MAi5,0,0) 


= DihLiSH 




'-'0 


G(MA29,0,1) 


= MxTx{ei)CihLiS 


-! /j 1 1 1 Vtj 


A(MBig,0,0) 
'-'0 


G(MBig,0,0) 


= Ta{ei)D2hLiS-R 


{ftl,l,l}R 


^(MSa,, 0,0) 


G(MB29,0,0) 


= T4ei)D2ahLiSIL 


{/ii,i,i}r 


A(MBg,0,0) 

<^0 


G(M Eg, 0,0)1 


^T4ei)C2xhLiS-R 


{/ii,i,i}r 




G{MEg, 0,0)2 


= Tz{ei)C2ahLlSIl 


— /12,1,i}r 




G{M Eg, 0,1)1 


^ M,mei)C2xhLiS 


{ftl,l,l}R 




G{MEg,0, 1)2 


= M^T,{ei)C2ahLiS 


— /12,1,i}r 



ii = L(ei + 62, 62 — ei) 
T^{el) = {S,G2«T(6i)} 



2ah — {E, C2z, C2a, G26, /, IC2Z, IC2a, IC2b} 
C2xh = {E, C2x, I, IC2x} 
C2ah ~ {E, C2a, I, IC2a} 

Mi = {E,tC2i} 



= {X)ij = SijXi. (5-5) 

These represent that the occupation numbers of electron for three atomic orbitals 
(■015 V'2, V'3) are Ai,A2,A3, where 

tpi = 4>lUii + ^2t^21 + '?!'3C^31, 
tp2 = 4>lUl2 + 02^^22 + (t>?>U32, 

tp3 = (f)lUi3 + ^2t^23 + (l>3U33. (5-6) 



We give three examples to show how each axial isotropy subgroup determines the 
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Table VI. Axial isotropy subgroup and its fixed point subspace of G\ 



R-rep 


axial isotropy subgroup 


fixed point subspace 


^(XAg,0,0) 
^0 


G(XAg,0,0)i = D2hL2SB. 
G{XAg, 0,0)2 = DihLaSR. 


■f/ll 1 1 rR 

{?ll,l,l + /12,1,i}r 


^(XSig,0,O) 



G(XBig, 0,0)1 = Tr(ei)C2zhi25'R 
G{XBig, 0, 0)2 = r^(ei + e2)D2ahL3SIi 


{ftl,l,l}R 

{/ll,l,l + /i2,l,l}R 


^(XBlg,0,l) 




G{XBig, 0, 1)1 = MyT^{ei)C2zhL2S 
G{XBig, 0, 1)2 = MyT^ei + e2)D2ahL3S 


■f/ll 1 ll-R 
L ^)^! ^ J R 

{/ll,l,l + /12,1,i}r 


A(XB2g,0,0) 




G(XB29, 0,0)1 ^T^{ei)C2yhL2SIl 

G(XZ^2„. 0.0)2 = r, (ei)T„(e2)C2„;,£:iS'R 


j/ii 1 iIr 

{/'i.i.i + /'2,i,;i}r 


^(Xi32g,0,l) 


G(XB2g,0, 1)1 = Mzn(ei)C2yhL2S 

G(XB2g,0, 1)2 = M,T^{ei)Ty{e2)C2ahL3S 


{/ii,1,i}r 

I'll, 1,1 + /12,1,i}r 


^(XB3g,0,0) 


G(XB3s,0,0)i =Tj,(ei)C2.fei2SR 

G(XB3g, 0,0)2 = T4e2)Ty{{ei)C2ahL3SIi 


{'ii.i.iIr 

{/ll,l,l + /12,1,i}r 




G{XB3g, 0, 1)1 = M,Ty{ei)C2.hL2S 
G(XS3g,0, 1)2 = M,T4e2)Ty{ei)C2ahL3S 


{/ii,i,i}r 

{/ll,l,l + /12,1,i}r 



£2 = i(2ei, 62), £3 = L(2ei, 262) 

ri(ei +62) = {E,C2iT{ei +62)}, i = z,x,y,a,b 

Ti{em) = {E,G2iT{ei)}, i = z,x,y,a,b, m = l,2 

Mi = {E, tC2i},i = z, X, y, a, b 

D 2ah — {E,C2z,C2a,C2b, I , IC2Z, IC2a, IC2b} 
C2xh = {E,C2x,I,IC2x} 
C2ah = {E,G2a,I,IC2a} 



canonical form of tlie HF Hamiltonian Hjn, occupied orbitals and their occupation 

numbers, and the type of the LOP of the state. 

Example 5.1. G{rAig, 0,0) state. 

In this case the isotropy subgroup G{Hjn) of Hjn is 

G{rAig,0,0)=D^hLoSIi. (5-7) 
Prom T(ei) and T(e2) G Lq invariance of we see that only is non-zero. 
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RP^ is invariant under C2z, C2x-, C2a £ -D4/1 and t R. Thus using p-19p we have 



C2z • R 



00 



C2x ■ R 



00 



C2a • R 



00 



^21 

pOO 


oOO 
-f>^12 
puu 
-"■22 
pOO 


-R'l^ ) 

-^23 , 

R's^ J 


= R"" 


pOO 

-R21 


pOO 
-"^22 
oOO 
-"32 


-RTs \ 

pOO 
-"23 

Rfs J 




pOO 
-f1^22 

RT2 

-Rf2 


pOO 
R2I 

RTi 
-RTi 


-RTs \ 
-R'l^ , 
^3°^ ) 


1 =i?oo 



t • = = (5.8) 



Thus we obtain 



= I a I , (5-9) 




where a and 6 are real numbers. From (I3-7P and SCF condition p-Sp we obtain 

X?? = N{Wu1iRTi + Wl212i?22 + VFi313i?33) 

= {U + 2U' - J)a + {2U' - J)b = a, 

= iV(H^212li??? + W2222RII + ^^2323^3°) 

= {U + 2U' - J)a + {2U' - J)b = a, 

= iV(T^313li??? + W^2^2Rll + H^3333^^!l) 

= 2{2U' - J)a + Ub = p. (5-10) 
Thus we obtain HF Hamiltonian Hm as follows: 

Hm = Hk + X]'L"(4fcsOlfcs + 4fcs«2fcs) + /?4fc^a3fcJ 
fc s 

= Hk + a/i(rAig, 0, 0)1,1,1 + fih{rA%,Q, 0)i,i,i. (5-11) 

Here the bases 0, 0)i,i,i(j = 1,2) are given in Table Hill This state corre- 

sponds to the normal paramagnetic state. 

Now we consider M point non magnetic states. All isotropy subgroups of M 
point non magnetic states contain Li = L[ei + 62,62 — 61) and S. Thus from S 
and Li invariance of R'^, only R^^^ and R}^ are non- zero. From ()3-26p . for m such 
that T{m) S Li, we have 

D{m) = D^\m) = D^^{m) = + 
D{m + ei) = D^\m + e^) = i:)^^(m + 6^) = - (5-12) 



Symmetry Classes of Spin and Orbital Ordered States in a t2g Hubbard Model 23 



where i = 1,2. Thus diagonahzing D{m) and D{m + e,), we obtain occupied 

atomic orbitals and their occupation numbers at sites m and m + ej. We consider 

the G(M52g,0,0) state. 

Example 5.2. G(MS2g,0,0) state. 

Since 

G(MS2g,0,0) =r,(ei)£)2a/,Li5R, (5-13) 

and R}^ are invariant under C2z , C2a £ D2ah and t & R. Thus using ()3-19p we 
have for / = 0, 1 



r< 7?'0 — pro nlO nlO \ _ 



C2a • = I I = 

t- R^'^ = {R'y = R'". (5-14) 

Then i?"' have forms 

i?™ = 1^ i?™ 41 I . (5-15) 

From C2xT{ei) invariance, we obtain 

/ ii?? \ 

V ; 

/ -iil? \ 

C2xTiei) ■R^^= R\\ -R\l \ = R^^ . (5-16) 

V -Rll 

Thus we obtain 

= I a I , = I c I , (5-17) 



^11 

--^31 


^12 
-^22 
--^32 


--^13 
-4°3 
^33 


^22 

^i°2 
p/0 
--^32 


--^31 


-i?^°3 
^33 


my ^ ^ 


10 




[f 


R[% 
R{\ 



") 

4°3 / 





where a, b and c are real numbers. From (j3-7|) and SCF condition (|3-8|) we obtain 

4°) 



X?0 = N{WnilRTl + W^1212i?22 + VFi313i?3°' 



= ([/ + 2?7' - J)a + {2U' - J)b = a, 

X^^ = N{W2l2lRu + W^2222i?22 + W^2323i?33) 

= {U + 2U' - J)a + {2U' - J)b = a, 
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xll = N(Wu3iRTi + W3232R'^2 + ^^3333^33) 



lOO 



00 



2{2U' - J)a + Ub = P, 



x\9, = N{Wn22Rl'l + Wi22lRl^2) 



10 



?10n 



™10 
^21 



(f + 2 J - U')c = 7, 

iV(VF221li?l^ + ^2112i?2?) 

= ( J' + 2 J - C/')c = 7. (5-18) 
Thus we obtain the HF Hamiltonian Hm as 
H^ = Hk + ^{a(,a\^,a^^s + 4fcs«2fcs) + /?4fc^a3fcs 

= + ah{rA\g,0, 0)1,1,1 + 0, 0)1,1,1 + -/h{MB2g,0, 0)i,i,i. (5-19) 

The fourth term of (I5T9P is the primary part which leads to the transition to the 
G(MS2g,0,0) state. From ([342]) we obtain the HF energy Euf as 

Erf = {HkY 

+ N{2{U + 2U' - J)a^ + Ub'^ + 4(2^7' - J)ab + 2(2J + J' - U')c^]. (5-20) 
For m such that T{m) G Li we obtain 



00 



R 



10 




D{m + ej) 

Diagonahzation of D{m) and D{m + e^) are written as 

U^D{m)U 

U^D{m + ej)U 



(i = l,2). 



(5-21) 




(5-22) 



where 



U 



1 

V2 

V 



1 

1 

V2 





\ 



1 / 



(5-23) 



From (j5-6p we obtain the occupied atomic orbitals and their occupation numbers for 
G(MB23,0,0) state as shown in Table EUl 
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Table VII. Occupied atomic orbitals and their occupation numbers for G{MB2g, 0,0) state 



site spin atomic orbital occupation number 



m 


n 


^1 = 




+( 


h) 


a + c 


m 


n 


1p2 = 




(^1 — ( 


h) 


a — c 


m 


n 


1p3 = 


03 






b 


m + Sj 


n 


V"! = 




+< 


h) 


a ~ c 


m + Sj 


n 


V'2 = 




(^1 — ( 




a + c 


m + Bj 


n 


V'3 = 


03 






b 



m : T{m) G Li, j = 1,2, a = ii?? = b = c = = i?£J? 



y 




occupation number a + c in the G{MB2g, 0, 0) state. 

In FigH] we show the pattern of orbital order by -^{'f'l + 4'2) and -^{(f'l ~ 
(f>2) with the occupation number a + c. Note that this pattern has symmetries of 
C2a,C2b,C2z,C2xT{ei) and C2xT{e2) G G{MB2gO,0). From ([3^ we see that this 
state has alternating quadrupole moments, for m such that T{m) G Li, as follows 

Qi2{m) = 4/2C, Qi2(m + ej) = -Ahc, j = 1, 2. (5-24) 

Finally we consider X point non magnetic state. As shown in Table IVTl in cases 
of G{X^,0,h')i, isotropy subgroups contain L2 = -L(2ei,e2) and S. Thus from 
L2 = L(2ei,e2) and S invariance of R''^, only R''^ {I = 0,2) are non-zero. From 
(j3-26p . for m such that T{m) G L2, we have 

D{m) = D^\m) = D^^{m) = + R^^ , 
D{m + ei) = D^\m + ei) = D^^{m + ei) = - R^^ . (5-25) 

In the cases of G(X7, 0,i/)2, isotropy subgroups contain Z3 = Z(2ei,2e2) and S. 
Thus from is = -L(2 ei,2e2) and S invariance of R^^, only {I = 0,1,2,3) are 
non-zero. From (j3-26p . for m such that T{m) G L2, we have 
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D{m + ei) = D^^m + ei) = D^^{m + ei) = - R^" - R^" + , 
D(m + 62) = r>"(m + 62) = £>^Hm + 62) = - Ji'" + fi'" - 
D{m + 61 + 62) = £)^^(m + 61 + 62) = D^'^im + 61 + 62) 

= R'^" + R''' -R^° -R^°. (5-26) 

We consider the G^XB-^g, 0,1)2 state which breaks time reversal symmetry as an 
example. 

Example 5.3. G{XB3, 0,1)2 state. 
Since 

G{XB3g, 0, 1)2 = M,T,{e2)Tyiei)C2ahL3S, (5-27) 
i?"^ {I = 0, 1, 2, 3) are invariant under G2a, C2xT{e2), G2yT{ei),tG2z- Thus we obtain 

/ a \ 
= a , R'' = 





\ 

i?20 = [ id \ , i?30 = I I , (5-28) 

-id ) 

where a, h, c and d are real numbers. Using (j3-8p we obtain i/^ as 

Hra = Hk + ^ ^{a{a\j^^a^hs + '^\ks"'2ks) + f^^lks'^'ihs 
k s 

•e/ t _ t 

*"Va2(fc+Q2)s"3fcs ^i{k+Q^)s°-2ks 

~ 4(fc+Q3)s"'3fcs + 4(fc+Q3)s"'lfcs)j' 

= Hk + ah{rA\g,0, 0)1,0,1 + /3/i(rA?g, 0, 0)1,0,1 + lh{MB2g, 0, 0)i,o,i 
5{h{XB3g,0, 1)1,0,1 + /i(X53g, 0, 1)2,0,1}, (5-29) 

where a, (3, 5 and 7 are determined by SCF conditions 

a = {U + 2U' - J)a + {2U' - J)b, 
(3 = 2{2U' - J)a + Ub, 
7 = (2 J + J' - [/')c, 

6 = {2J - J' - U')d. (5-30) 

The fifth term of (j5-29p is the primary part which leads to the transition to the 
G(Xi?3g, 0, 1)2 state. The fourth term of (j5-29p is the secondary part induced by 
the transition. 

The HF energy is written as 

Shf = (^k)' + y {2(f^ + 2f/' - J)a^ + + 4(2^7' - J)ab} 
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+ y{2(2J + J' - U')c'^ + 4(2J -U' - J')d^}. (5-31) 

From (15-26P we obtain, for m such that T{m) G L3, 

I a c —id \ 
D{m) = c a id \ , D{m + 62) = 



a 


c 


—id 


c 


a 


id 


id 


—id 


b 


a 


—c 


—id 


—c 


a 


—id 


id 


id 


b 




D{m + ei) — —c a —id , D{m + ei + 62) = c a —id . (5-32) 



Diagonahzing D{m), D{m + ei), D{m + 62), D{m + ei + 62), we obtain occupied 
Table VIII. Occupied atomic orbitals and their occupation numbers for G{XBzg, 0, 1)2 state 



site 




spin 


atomic orbital 


occupation number 


m 
m 
m 




Ti 

n 
n 


^1 = 

l/'2 = 
lAs = 


^sgn(d)ii(;ii - ^sgn(d)ii(?!)2 + iwcps 
-^sgn{d)w(j)i + ^sgn(d)w<^2 + iu03 


Ai 
A2 
A3 


m + ei 
m + ei 
m + ei 




n 
n 
n 


ii2 = 

tp3 = 


-^sgn(d)u0i - ^sgn(ii)M02 - iw(/>3 
-^sgn(d)u;(^i - ^sgn(d)w(?!)2 + m^s 
^(<^i - <^2) A3 


Ai 
A2 


m + 62 
m + 62 
m + 62 




n 
n 
n 


= 

'4>2 = 
1p3 = 


-^sgn{d)u(j>i + -^sgn{d)u(t>2 - iw(j)3 
^sgn(d)iu<^i + ^sgn(d)«;02 + mc^ia 
-^(4>i-4>2) 


Ai 
A2 
A3 


m + 61 
m + 61 
m + 61 


+ 62 
+ 62 
+ 62 


n 
n 
n 


Ipl = 
^3 = 


-^sgn(d)u0i + ^sgn(d)M02 + iw<j)-i 
^sgn(d)TO(^i - ^sgn(d)ui<;!>2 +iu(j)3 


Ai 
A2 
A3 



r(m) e L3 = (261,262) 

„ _ rOO _ pOO > _ pOO „ _ plO _ plO -I _ p20 _ p20 _ p30 _ p30 
a — itii — -rt22, O — ^331 C — -K12 — -rt21, — ^23 — "^32 — —^13 — -"-31 



w= -4(1+ ; )Kw = ^(1 



, a + b-c+y/(a-b-cp+8d^ . a + 6- c- ^/("i-ii-c)^ +8d2 

Ai = ^ — 2 , A2 = — 2 , A3 = a + c 

atomic orbitals and their occupation numbers at the sites m, m + ei, m + 62, m + 
ei + 62 for the G(Xi?3g, 0, 1)2 state as shown in Table IVlnl Note that occupied 
atomic orbitals have complex coefficients. 

From (j3-43p . we see that this state has the following orbital angular momen- 
tum(the primary LOP) and quadrupole moment(the secondary LOP). For m such 
that T[m) € L3, we obtain 



li{m) = 4d, him) = id, 
li{m + ei) = —M, him + ei) = 4(i, 
/i(m + e2) = 4d, h{m + = -M, 
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©1+62 






Fig. 2. The ordering pattern for quadrupole moments: Q12 and orbital angular momenta: 
I = {li,h) in the G{XB3gO,l)2 state. The ovals and arrows represent Q12 and orbital an- 
gular momentum I, respectively. 



The ordering pattern of Q12 and / = (/i, I2) for the G{XB3g,0, 1)3 state is shown in 
Fig. O Note that the onset of quadrupole moment Q12 is induced by the transition 
to the G{XB3g,0, 1)2 with the orbital angular momentum. 

By a similar manner to the case of G{XB^g, 0,1)2, we can see that all states 
having isotropy subgroups such as G{A,0,1) or G{A,0,l)j (j = 1,2), which break 
the time reversal symmetry, have complex occupied orbitals and orbital angular 
momentum as LOP. 

§6. Symmetry classes of magnetic orbital ordered states 

In order to list magnetic orbital ordered states with ordering vectors Q; {I = 
0, 1, 2, 3), bifurcating from the normal state, we present the axial isotropy subgroups 



of the R-reps Gq '^'^''^ ,Gq and Gq , The axial isotropy subgroups of R- 

reps gJ,^^'^'''^ and gJ,^^^'^'''^ are hsted in Table US Those of g'^^^'^'"^ are listed in 



Spin magnetic states are classified into two groups: collinear and non-collinear 
magnetic states. The collinear magnetic state has the isotropy subgroup containing 
the subgroup ^(63): the spin rotation around z axis. The isotropy subgroup of the 
non-collinear magnetic state does not contain A{e-j). 

Non-collinear magnetic states are derived from the R-reps: G^^^^'^''^\ G^^^^''^'^\ 
q{xAi„i,i)^ ^{XBjg,i,i) ^{XBjg,ifl) ^^^^^ A = r,M and j = 1,2,3. The corre- 
sponding non-collinear magnetic states are G{rEg, 1, 1)3, G{rEg, 1, 0)3, 



G{MEg, 1, 1)3, G{MEg, 1, 0)3, G(X^3, 1, 1)3, G{XBig, 1, 1)3 and G{XBig, 1, 0)3 (i = 



li{m + ei + 62) 
<3i2(m) 
Qi2(m + e2) 



-4d, him + ei + e2) = -M 
AI2C, Qi2{'m + ei) = -4/2C, 
-4/2C, (5i2(m -h ei + 62) = 4/2C 



(5-33) 



Table m 
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Tal)k> [X. 


Axial isotropy 


subgroup and its Fixed point subspaco of fr',',^' 


and G[^'-'"-"' 




axial isotropy subgroup 




^(^^413, 1,1) 


G(r^i„i,i) = 


= M(62)£>4hioA(63) 


{hi,i,i}B. 




G(r^2s,i,o) : 


= D4h(£;,M2a:)ioA(e3)R 


{/ii,3,i}r 


^(rsig,i,i) 


G(ms,i,i) : 


= M{e2)D4h{u2a:,E)LoA{e3) 


{ftl,3,l}R 


/,(rB2g,l,l) 


G(rB29,i,i) = 


= M (62 )L>4h('U2i:,U2a;) £0^(63) 


{/ii,3,i}r 




G(riJ„ l,l)l : 


= M(62)Z32h(u2:., £;)LoA(63) 


{/ii,3,i}r 






= M{e2)D2ah(u2x, E)LoA{e3) 


{hl,3,l + ft2,3,l}R 




G{rEg, 1,1)3 ■ 


= M{e3)D4h{ut^,U2x)LQ 


{hl,l,l + /12,2,i}r 


^(rB„ ,1,0) 


G{rEg, 1,0)1 ■ 


= D2h{u2-^,E)LoA{e3)Ti 


{/ii,3,i}r 




G(rEg, 1,0)2 ^ 


= D2ah{'U.2a:, E)LoA{e3)R 


{ftl,3,l + ft2,3,l}R 




G(rEg, i,o)i ■- 


= D4h{u4^,U2x)LoR 


{ftl,3,l + fe,2,l}R 


^ f A/f A , 1 1 \ 


G{MAlg, 1,1) 


= M(e2)r"(6l)£>4hilA(63) 


{/ii,i,i}r 


P=,(MA2g,l,0) 


G{MA2g,l,0) 


= T,(6i)T"(6i)£>4/.(£^,W2a=)ilA(e3)R 


{/Ii,3,i}r 


^(MSlg,l,l) 

'-'0 


G{M Big, 1,1) 


= M(e2)T='(ei)Ta(ei)i34h(M2., A(e3) 


{/11,3,i}r 


P=,(MS25,1,1) 


G(MS2s,l,l) 


= M(e2)T'^(ei)T^(ei)r>4h(w2x, W2x)ii A(e3) 


{/11,3,i}r 


A(MiSg, 1,1) 


GiMEg, 1, 1)1 


= M{e2)T''{ei)Ty{ei)D2h{u2a:,E)LiA{e3) 


J /l 1 01 

1 1)3, 1 jrt 




G{MEg, 1, 1)2 


= M{e2)T^ei)n{ei)D2ah{u2.,E)LiA{e3) 


{'11,3,1 + '12,3,i}r 




G{MEg, 1, 1)3 


= M{e3)T:,{ei)Dih{ut^,U2x)Li 


{/ll,l,l + /12,2,i}r 


A(MBg, 1,0) 


G(M£g, 1,0)1 


= T''iei)Ty{ei)D2h{u2., E)LiA{e3)Il 


{/ii,3,i}r 




G{M Eg, 1,0)2 


= T''{ei)Ty(ei)D2ah(n2:., E)LiA{e3)K 


{/ll,3,l + /12,3,i}r 




G{MEg, 1,0)3 


= r"(6l)r,(6l)£>4h(wi,M2x)ilR 


{/ll,3,l + /12,2,i}r 


io = i(ei, 62 


;), Ll = i(6i + 62, 62 - ei) 




M2i = u{ei, tt) 


G S, M{ei) = 


{E,tU2i} 





A(ei) = {^(63,6') \ 0<e<2Tv} 

D4h{a,l3) = {(£;,G4+a,G2.a',G4-a-) + G2./3(S,G4+a,G2za",G4-Q-)} x Ci 
D2h(a,l3) = {(£;,G2.a) + G2./3(£;,G2.a)} x Ci 
D2ah{a, 13) = {{E, C2.a) + C2al3{E, C2.a)} x Ci 
T'iem) = {E,U2iT{em)} 

Tr{em) = {E,C2rT{era)} 
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Table X. Axial isotropy subgroup and its Fixed point subspace of Gq^'''^''^'' 
R-rep axial isotropy subgroup Fixed point subspace 



^(XAlg.l.l) 

'-'0 


G{XArg 


1, 1)1 = M{e2)T''iei)D2hL2A{e3) 


{^1,3,i}r 






G{XAi, 


1, 1)2 = M(e2)r"(ei + e2)D4hL3A{e3) 


{/ll,3,l + ^12, 3, 


i}r 




G{XAlg 


1, 1)3 = M(e3)r"(e2)T^(ei)£)4h(u2a, E)L3 


{hl,l,l + h2,2, 


i}r 




G{XBlg 


1, 1)1 = M{e2)T''{ei)D2h{E,U2^)L2A(e3) 


{^1,3,i}r 






G{XBi, 


1, 1)2 = M(e2)r"(ei + e2)D4h{u2^,U2^)L3A{e3) 


{/ll,3,l + ^12, 3, 


i}r 




G{XBrg 


1, 1)3 = M(e3)r"(e2)r''(ei)£)44^26,^2.)£3 


{/ll,l,l + ^2,2, 


i}r 


"-'0 


G{XB2g 


1, 1)1 = M(e2)T==(ei)D2h(M2:.,M2.)£2A(e3) 


{^1,3,i}r 






G{XB2g 


1, 1)2 = M{e2)T^{ei)Ty{e2)D2ah{u2a^,E)L3A{e3) 


{/ll,3,l + ^12, 3, 


i}r 




G{XB2g 


1, 1)3 = M(e3)r,(ei)Ty(e2)i34h(Mt,^t2.)i3 


{/ll,l,l + ^2,2, 


i}r 


"-'0 


G{XB3g 


1, 1)1 = M{e2)Ty{ei)D2h{u2^, E)L2A{e3) 


{^1,3,i}r 






G{XB3g 


1, 1)2 = M(e2)r,(e2)T„(ei)i)2ah(u2.,S)i3^(e3) 


{hl,3,l + h2,3, 


i}r 




G{XB3g 


1, 1)3 = M(e3)r,(e2)T^(ei)i)4h«,,^t2.)i3 


{hl,l,l + h2,2, 


i}r 




G{XBlg 


1, 0)1 = T^{ei)D2h{E, U2^)L2A{e3)Yi 


{ftl,3,l}R 






G{XBrg 


1,0)2 =T"(ei +e2)L»4h(u2^,U2:.)£3A(e3)R 


{hl,3,l + h2,3, 


i}r 




G{XBlg 


1,0)3 =T"(e2)T^(ei)i)4h(«26,W2.)i3R 


{/ll,l,l + /l2,2, 


i}r 


^(XS2g.l,0) 

"-'0 


G{XB2g 


1,0)1 =T'"(ei)£>2h(M2.,^*2.)i2A(e3)R 


{'ii.s.iIr 






G{XB2g 


1, 0)2 = T4ei)Ty{e2)D2ah{ii2^,E)L3A{e3)'R 


{^11,3,1 + ^12, 3, 


i}r 




G{XB2g 


1,0)3 ^ T4ei)Ty{e2)D4h{ut,,U2:c)L3R 


{^1,1,1 + ^12, 2, 


i}r 


"-'0 


G{XB3g 


1,0)1 =T'j,(ei)D2h(M2.,£)-L2A(e3)R 


{^1,3,i}r 






G{XB3g 


1, 0)2 = T,(e2)rj,(ei)£)2ah(u2,, £)£3A(e3)R 


{^1,3,1 + /l2,3, 


i}r 




G{XB2g 


1,0)3 = T4e2)Ty{ei)D4h{ut,,U2x)L3R 


{/ll,l,l + /«2,2, 


i}r 


L2 = L{2ei, 


62), L3 = 


£(261,262) 






U2i = u{ei,-n 


),U2a — u{ei + 62, Tt), W2i, = u( — 61 + 62, Tt) G S 






M{ei) = {E,tu2^}, T\ 


6^) = {s,u2.r(e„)}, T,(6„) = {£;,C2»r(6„)} 






D2h{a,l3) = 


{{E,C2,a 


) + C2./3(S,C2za)} X Ci 






D2ah{a, P) = 


= {{E, C2.Q) + C2af3{E, C2.a)} X Ci 






D4h{a,/3) = 


{{E,C+a 


,C2zO?,C4^a ^) + C2xP{E,C4^a,C2za'^ ,C4^a ^)} x Ci 





1, 2, 3). All states in Table HXl and 1X1 except these eleven states are collinear magnetic 
states. We consider collinear and non-collinear magnetic states seperately. 
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6.1. Collinear magnetic state 

All axial isotropy subgroups of collinear magnetic states contain ^4(63), then 
j^ii ^ ^ Q / = 0, 1, 2, 3. Thus in these cases we obtain from (l3^ 

D'^^im) = D^^m) = 0. (6-1) 

We consider two examples. 

Example 6.1. G{rAig, 1, 1) state. 

In this case, the isotropy subgroup G{Hm) of Hm is 

G{rAig, 1, 1) = M{e2)D^hLoA{es). (6-2) 

From Lq invariance of ii'^ we can see that only RP^ and RP^ are non-zero. 
From C2X, C2Z, C2a{^ -D4/1) invariance of R^^ and R^^, we obtain 



a 





M 




( - 










a 










c 










b 1 




I 





d / 



= I a I , i?"^ = I c I , (6-3) 

\006/ \00: 

where a, 6, c and d are real numbers. From (j3-26p we obtain 

D^i(m) = i?°°-i?°3 = I a-c 1. (64) 

\ 6-d / 

In Table IXH we list the occupied atomic orbitals and their occupation numbers for 
G{rAig, 1, 1) state. From (I3-43P the spin density at the site m is given by 

Table XI. Occupied atomic orbitals and their occupation numbers for G{rA\g, 1, 1) state 




site spin atomic orbital occupation immber 



m 


T 


i>i = 


01 


a + c 


m 


T 


i>2 = 


(t>2 


a + c 


m 


T 


i>-i = 




b + d 


m 


i 


ipi = 




a — c 


m 


i 


i>2 = 


02 


a — c 


m 


i 


i>-i = 


03 


b-d 


a = 


pOO 
itll 




pOO „ _ p03 
^33 1 C — itii 


_ p03 J _ p03 






s^{m) = 




2c + d 



This state corresponds to the usual ferromagnetic state without orbital order. 



(6-5) 
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Next we consider M point collinear magnetic state. From Li invariance of i?''^, 
only r!'^{1 = 1, 0, a = 0, 3) are non-zero. 
Example 6.2. G{MBig, 1,1) state. 
In this case, the isotropy subgroup G{Hm) of Hm is 



G{MBig, 1, 1) = M{e2)T^{ei)Ta{ei)D4h{u2., E)LiA{e3). (6-6) 



From C2x, C22; C2aT{ei) € G{MBig, 1, 1) invariance, we can see that only RP^ and 



R^"^ are non-zero and 

= I a I , i?^^ = I -c I , (6-7) 
\006/ \000/ 

where a, b and c are real numbers. From ()3-26p . we obtain for m such that T{m) £ Li 

/ a + c \ 
D^\m)=\ a-c , i:>^^(m) 







D^^m + ej)=\ a -he , i:)-^-^(m) = a-c .(6 



In Table IXIII we list the occupied atomic orbitals and their occupation numbers for 
G{MBig, 1,1) state. From (j3-43p . we see that this state has the spin quadrupole 

Table XII. Occupied atomic orbitals and their occupation numbers for G{MBig, 1, 1) state 



site 


spm 


atomic orbital 


occupation number 


m 


T 


in = 1 


h 


a + c 


m 


T 


1p2 = ' 




a — c 


m 


T 


l/)3 = , 




b 


m 


i 


Ipl = ' 




a ~ c 


m 


i 


V'2 = ' 




a + c 


m 


i. 


^3 = ' 




b 


m + gj 


T 






a — c 


m + Bj 


T 


V'2 = ' 


^2 


a + c 


m + Bj 


T 


ips = . 


i>3 


b 


m + gj 


i 


Ipl = i 


h 


a + c 


m + gj 


i 


1p2 = ' 


^2 


a — c 


m + gj 


i 


ipz = ' 


^3 


b 



j— 1,2, a — Till — ^22, b — i?33, c — Rli — —R22 
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moment Q2 {tti), for m such that T{m) G Li, as follows: 



Qf{m) = 2c/i,Qf (m + e 



23/ 



-2cli,i = l,2. 



From (13-43P we obtain 

Ql{m) 

s^{m) 



R22) 



13 



c = 0, 



0. 



(6-9) 



(6-10) 



Thus G(Mi?ig,l,l) state, which has the spin quadrupole moment ■, is not the 
coexisting state of quadrupole moment Q'^iyn) and spin density s^(m). 

6.2. Non-coUinear magnetic state 

First we consider M point non-collinear magnetic state. 
Example 6.3. G{MEg, 1,1)3 state. 
In this case, the isotropy subgroup G{Hm) of Hm is 

G{MEg, 1,1)3 = Mie3)T,iei)Dihiut,U2,)Li. 



From T(ei + 62) G Li invariance of i?''^, we can see that R 



,2A 



R 



3A 



(6-11) 
0. From 

tu2z G M(e3) invariance, we see that and = 0,1) are real matrices 

and R^^ and /^^'^ are pure imaginary matrices. From C2zT{ei) and C2zU2z,C2xU2x, 
C2aU2a £ D ^^{u'l^ , U2x) invariancc we can see that only R^^,R^'^,R^^ and R^^ are 
non-zero and have following forms. 

a 



R' 



lOO 








a 




i03 



-IC 




R 



11 




d 
d 



R 



12 








-d 









(6-12) 



where a, b, c and d are real numbers, and i is imaginary unit. From (j3-26p and 
p-29p . we obtain, for m such that T{m) G Li, 



T>{m) 



■D{m + Bj) 



I 



a 

-ic 




-id 



a 
—ic 




y id 



ic 
a 



d 

ic 
a 








b 

-id 
d 




b 

id 
-d 







id 
a 
ic 







d 

-ic 
a 




id \ 
d 



b J 






—id 
a 
ic 







-d 
-ic 

a 





-id \ 
-d 







b 



(i = i,2). 



(6-13) 
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Diagonalizing T>, we obtain occupied general spin orbitals and their occupation num- 
bers as in Table IXIIII From (|3-43p we see that this state has following local order 

Table XIII. Occupied general spin orbitals and their occupation numbers for G{MEg, 1, 1)3 state 



site 


general spin orbital 




occupation number 


7n 


= 


^^Z<^l+02) IT) 






fii 


l/'2 = 


^w(-z<^i + 02) li) + 


i3 It) 




ffi 


1p3 = 


2^«;H0i-02) |T)+«<? 


k 1 1 \ 

^3 14-/ 


Xo 


'lit 


^4 = 


^l^(i01-02) |i)+U03 


1 1 / 


A2 


m 


V'5 = 


^(01+i02) IT) 




A3 


m 


1p6 = 


TsWi-iV'a) U) 




A3 


m + Bj 


V"! = 


^u(z0i + 02) !T) -i^-^a 


li) 


Ai 


m + Sj 




^w(-i0i + 02) U)-m;<7 


^3 IT) 


Ai 


m + Bj 


ips = 


^^H0i-02) IT)-"-? 


'3 li) 


A2 


m + ej 


■04 — 


2S^«;(i0i - 02) U)-^03 


IT) 


A2 


m + gj 


V'S = 


75Wi+i<^2) IT) 




A3 


m + gj 


V'6 = 






A3 



T{m) e Li 

7 = 1,2, u=^(l+ , )i, w = ^(l- , )i 

„ — pOO _ pOO L _ nOO „• _ p03 _ p03 , _ pll _ pll _ ol2 _ pl2 
a — till — -K22, O — _K33, ZC — tii2 — —1X21,0. — U23 — -K32 — —^13 — —^31 

Ai = ^^^^5 , A2 = ^^^^5 , A3 = a - c 



parameters. For m such that T{m) £ Li and j = 1, 2, we obtain 

Qlaim) = 2l2d, Qlsim + ej) = -Ihd, 

Ql^{m) = -2l2d, Qliim + ej) = 2/2^, 

/|(m) = /|(m + Bj) = 2c. (6-14) 

Note that the onset of ferro spin obital angular momentum:^! is induced by the transi- 
tion to the G{MEg, 1, 1)3 state having non-collinear spin quadrupole moment:{Q235 Qsi}- 
From ()3-43p we can see that for all site m 

s^(m) = s^(m) = s^(m) = 0, 
him) = 0, 

Q23(m) = Q3i(m) = 0. (6-15) 

The important point to note is that the existence of order parameters /^(m), (523 (^) 
and Q'^i{m) does not mean the coexistence of spin density s^(m) and orbital an- 
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gular momentum /3(m) nor spin densities s^{m),s'^{m) and quadrupole moment 
Q23{m),Q3i{m). 

Next we consider X point non-collinear magnetic state. We consider two exam- 
ples. 

Example 6.4. G{XAg,l,l)s state. 

The isotropy subgroup of the G{XAg, 1, 1)3 state is 



G{XAg, 1, 1)3 = M{es)T^{e2)Ty{ei)D^h{u2a,E)Ls. (6-16) 

From U2xT{e2) (G T^'(e2)) and U2yT{ei) (g T^(ei)) invariance of B}^, we see that 
only R^^,B}^,R^'^ and E?"^ are non-zero. From D2h-, C^zU2a (G Dih{u2a,E)) and 
M(e3) invariance of HI'^ we see that only RP^,R^^ and R^"^ are non-zero and have 
following forms, 

/ a \ /cOO\ / d \ 

i^oo = a \ , R^^ = i d \ , R^^ = i c , (6-17) 
\0 b J \0 f J \oo// 



where a, b, c, d and / are real numbers. From (13-260 we obtain 



D^^{m) = R 



.00 



D^^ (m) = e 



iQi Tn r>21 



R^' + ie 



(6-18) 



Then from (I3-29|) we obtain, for m such that T{m) G ^3, 



■D{Tn) 



a 



c + id 






a 



d + ic 







b 






c — id 


a 







d — ic 


a 







f-^f 


b 



■Dim + ei) = 



a -(c + id) 

a -(d + ic) 

& + 

-(c-id) a 

-(d-ic) a 

-{.f-if) 0b 



■D{m + 62) = 



a 








c + id 











a 








d + ic 











b 








f + ^f 


— id 








a 











d — ic 








a 






\0 f-if b J 
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Table XIV. Occupied general spin orbitals and their occupation numbers for G{XAg, 1, 1)3 state 



sit6 


general spin orbital 


occupciTjiori iiuriiuci 


m 


tpi = 


^(IT>+e'^ |i» 


Ai 


m 


1p2 = 


±|(|t)+e'(f-^) li)) 


Ai 


m 


V'3 = 


^(|T)+e'(-+^) U» 


A2 


m 


tp4 = 


^(|T) + e'(*-^') li)) 


A2 


fii 


i>5 = 


^(IT) + e^ li)) 


A3 


m 


i>6 = 


l|(||)+e'(-+i) li)) 


A4 


Tfl -\- Cl 


tpl = 


^(|T)+e'(--^) li)) 


Ai 


171 -\- ei 


1p2 = 


^(||)+e'(t+^) li)) 


Ai 


m + ei 


ti^ 


^(IT) + c-'^ li)) 


A2 


m + ei 


^4 = 


^(IT)+e'(-*+''' li)) 


A2 


TTi + ei 


V'5 = 


i|(||)+e«U-f) ID) 


A3 


m + ei 


V'6 = 


^(IT)+e'(-^> li)) 


A4 


771 + 62 


Ipl = 


^(IT)+e'(-^> li)) 


Ai 


771 + 62 


V'2 = 


i|(||)+e*(-f+7) ID) 


Ai 


77X + 62 


1p3 = 


i|(||)+e«(-7) ID) 


A2 


771+62 


Tpi = 


^(|T)+e'(t+^) li)) 


A2 


77i + 62 


V'5 = 


^(IT)+e''-^' li)) 


A3 


771 -|- 62 


V^e = 


^(||)+e»(-f) ID) 


A4 


771 + 61+62 


V*! = 


4|(|D+e'(-+^) ID) 


Ai 


771 + 61+62 


V'2 = 


^(|D+e*(-f-7) ID) 


Ai 


771 + 61+62 


V's = 


^(IT)+e^^ li)) 


A2 


771 + 61+62 


V'4 = 


^(IT)+e'(t-^' ID) 


A2 


771 + 61+62 


^5 = 


^(ID+e^U+f) ID) 


A3 


771 + 61+62 


V'6 = 


^(IT)+e''^' li)) 


A4 



r(77i) e L3, Ai = a + + d^ A2 = a - ^/cF+lP, \a = b + y/2f, A4 = & - \/2/ 
c + id 00 



Vc^ + 



pOO I _ dOO „ _ p21 _ p32 J _ p32 _ p21 r _ p21 _ p32 
H22, — C — Hii — -K22, " — rtii — -K22, / — -K33 — -K33 
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X>(m + 61+62 



a 



-(c + id) 







a 




-(d + ic) 







b 



-if + if) 



-(c — id) 


a 







-(d — ic) 


a 




\ 



-if^if) 


b 

(6-19) 



Diagonalizing J?, we obtain occupied general spin orbitals and their occupation 
numbers as in Table IXTVl 

In FiglS] we show the spin density pattern of s{m) = (s^(m), s^(m)) for four 
sites in the unit cell. Note that this spin density pattern has G{XAg, 1, 1)3 symmetry. 



Fig. 3. Spin density s(m) = (s^(m), s'^(m)) of the G(Xyl3, 1, 1)3 state. \s(m,)\ = c + d + f . 



Example 6.5. G{XBsg, 1,1)3 state. 

The isotropy subgroup of the G{XB3g, 1, 1)3 state is 



G{XB3g, 1,1)3 = M{e3)T^{e2)Ty{ei)Dih{ut,,U2x)L 



3- 



(6-20) 



From G{XB3g, 1, 1)3 invariance, we obtain non-zero R^\R^^ and as follows 



00 



R 



.21 



a 

a 

6 



d 

d 



R 



13 



ic 
-ic 




R 



,32 





-d 



-d 





(6-21) 



(6-22) 



where a, b, c and d are real numbers, and i is imaginary unit. From ()3-26p and (j3-29p 
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we obtain, for m such that T{m) £ 



X»(m) 



I>(m + ei) 



X>(m + e2) 



X>(m + ei +62) 



/ 


0. 


ic 


n 
u 








id 


\ 








— zc 


CL 


n 
u 








d 










n 
u 


n 
u 





id 


d 



















—id 


a 


-ic 
















d 


ic 


a 











v 


—id 


d 











b 


/ 






/ 


G 


— tc 


n 
u 










id 


\ 




IC 


fl 


u 












d 






U 


n 
u 


U 




id 


-d 















—id 


a 




ic 















-d 


—ic 


a 







v 


—id 


-d 













6 


/ 


/ 


a 


— ic 


n 
U 








-id \ 








ic 


a 


n 








d 










U 


U 


I, 



-id 


d 





















a 


ic 



















d - 


-ic 


a 











V 


id 


d 











b 


) 






/ 


a 


ic 













-id 


\ 




—ic 


a 













— 


d 












b 


—id 




-d 

















id 


a 




-ic 

















-d 


ic 




a 









V 


id 


-d 













b 




/ 



(6-23) 



Diagonahzing "D, we obtain the occupied general spin orbital and their occupation 
numbers as shown in Table IXV[ From (|3-43p we see that this state has spin orbital 
angular momenta l^{'m) and spin quadrupole momenta Q\'^{rn),Q\'^{rn) as follows: 



Ql 



23 



m 



Q\2,{m + 61 + 62 
ll{m 
/f(m + 6i 



Q 



m 



-2/2^, 



Q|i(m + 62) = 



-2/2d, 
-2hd, 



2hd, 
-2hd, 
2hd, 

-2l2d, Qliim + 61 + 62) = 2/2^, 
ll{m + 61 + 62) = 2c, 
/|(m + 62) = -2c, 



(6-24) 



for m such taht T{m) G -L3. Here we note that the secondary LOP:/! is induced by 
the appearance of the primary LOP:{Q23) Qii}- Prom (j3-43p we can see that for all 
site m 



s^{m) = s'^{m) = s'^{m) = 0, 



Symmetry Classes of Spin and Orbital Ordered States in a t2g Hubbard Model 39 

Table XV. Occupied general spin orbitals and their occupation numbers for G(XBzg, 1, 1)3 state 
site general spin orbital occupation number 



TTt 






,-sgn(d) / , ■ , X ||\ 1 , 


1 1 \ 

14-/ 


^1 


TTl 








•^•i 1 1 / 




TTt 








>3|i> 


A2 


TTt 








IT) 


Xo 
A2 












^3 


TTl 










A3 




ipi 




iM^w(01+j02) |T>+W03 


li) 




TTt ~\~ 61 


1p2 






is IT) 


Al 


TTt ~\~ Gi 


'03 






/\2 


TTt ~\~ Ci 








IT) 


Ao 
/\2 


TTt ~\~ Gi 


V's 








A,> 
A3 


m + Gi 






^(</.i+i«|i) 




A3 


TTt ~\~ G2 


l/jl 




+ i02) |T>+w<? 


isli) 


Xi 
^1 


TTt ~\~ G2 


1p2 




- i02) li) + W03 


IT) 






V'3 






li) 


Ao 

/\2 


TTt -\- G2 


1^4 






is IT) 


Xo 
A2 


TTt ~\~ G2 






75(<^i-^'^2) |T> 




Xo 

^3 


m + G2 










A3 


m + Gl +G2 






-e-s^u{cj>,-i<f>2) \})+W4 


isli) 


Ai 


m + Gl +G2 


V'2 




i?S^u{<i>i + icj>2) |i>+W03 


IT) 


Ai 


m + Gl +G2 


V'3 




- i</.2) |T)+«03 


li) 


A2 


m + Gl +G2 


l/)4 




-i2S^«,(01 + i02) |i)+«<? 


is IT) 


A2 


m + Gi + G2 


^5 




75('^i+*'^2) IT) 




A3 


m + Gl + G2 


V'e 








A3 



m : T{m) G Ls 

a = Rii = R22, b = i?33, ic = RI2 = —R21, d = R23 = R32 ~ —R13 ~ —Rl 

, a+b+c+\/(a-b+c)'^+ScP . a+b+c-\/(a-b+c)^+8<P . 

Ai = — 2 1 A2 = — 2 ' A3 = a — C 

— _l i^b±^^^\^ _ 1/1 a-b+c \^ 
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h{m) = 0, 

Q23(m) = Q3i(m) = 0. (6-25) 

Thus the G{XB-ig, 1, 1)3 state is not a coexistent state of spin density and quadrupole 
moment. 

Finally we note that many states in Table IIXI and |X] have spin orbital angular 
momentum and/or spin quadrupole moment, however they are not coexistent states 
of {spin density}, and {orbital angular momentum or quadrupole moment} except 
G{XAg, 1, 1)1 state. 

§7. Some numerical results 

In this section we present some calculated results of HF equations for the states 
listed in Table |Vl Table IVH Table HX] and Table Kl We solve the Hartree-Fock 
equation of a state (characterized by an isotropy subgroup Gj) self-consistently by 
starting initial values of R!"^ with Gj symmetry. From the SCF condition (j3-8p we 
obtain the initial values of a;'^. After diagonalizing Hm of (j3T0p with these a;''^, 
we obtain new The obtained i?'^ are substituted into (j3-8p to compute new 

x^^. We use them as inputs to repeat the above process until the relative errors in 
R^^ between successive iterations are less than the desired accuracy. After obtaining 
converged R}^ , we obtain the HF energy -Ehf from (|3-12p . 

We consider the cases of parameter sets listed in Table IXVII Since it is known 
that Ca2-a;Sr2.Ru04 has four 4d electrons in the t2g orbitals,^^'^-* we use number of 
electrons per site = 4.0 for all parameter sets. The eighth column denotes the 
most stable state among all states listed in listed in Table El Table IVH Table HX] and 
Table B 



Table XVI. Parameter sets for model calculations 



N.O 


5 


tl 


t2 




U 


j = j' 


most stable state 


(1) 


0.0 


1.0 


0.0 


1.0 


9.0 


2.25 


G(rAig,i,i) 


(2) 


0.4 


1.0 


0.0 


1.0 


9.0 


0.4 


G{MBlg, 1,1) 


(3) 


0.0 


1.0 


0.0 


1.0 


9.0 


0.7 


G{M Eg, 1, 1)3, G{X Bag, 1,1)3 


(4) 


-0.14 


1.0 


0.8 


0.8 


8.0 


1.0 


G{XAg,l,l)3 


number of electrons 


per 


site = 


4- 


4.0, U' = 


U-2J 



7.1. Parameter set (1) 

In the parameter set (1), where J has rather large values, the G{rAig, 1, 1) state 
is most stable. The calculated values of a, 6, c and d in (j6-3p are 

a = 0.6679, 5 = 0.6643, c = 0.3321, d = 0.3357. (7-1) 

Thus we obtain o + c ~ 1.0, a — c~ + d 1.0 and b — d i. Then from Table 
IXH the occupation numbers for ||) and (pi \[) are 



n{(t>i \]))=n{<P2 |T» =n(</.3 |T» = 1, 
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n{cPi\l))=nicP2\i))=n{cPs\i))^-. 



(7-2) 



7.2. Parameter set (2) 

In the parameter set (2), the G{MBig, 1, 1) state is most stable. The calculated 
values of a, 6 and c are 



Thus we have 



a = 0.5, h = l, c = -0.473. 



a + c = 0.027 « 0.0, a - c = 0.974 1.0, 6 = 1.0. 



(7-3) 
(7-4) 



From Table [Xm we obtain a qualitative pattern of the electron occupations as shown 
in FiglH 



f f tf 

^1 ^2 ^3 



01 02 03 



ei 

Fig. 4. Schematic pattern of the electron occupation in G{M Big, 1,1) state. | (J.) indicates one 
electron with up(down) spin. 

7.3. Parameter set (3) 

In the parameter set (3), the G{MEg, 1, 1)3 and the G{XB^g, 1, 1)3 states have 
the same energy and are most stable. The calculated values of a, h, c and d in (j6T2p 
and ([6^ are 

a = 0.7113, h = 0.5775, c = -0.2775, d = 0.3240. (7-5) 
Then the occupation numbers for general spin orbitals ipj {j = 1, ... ,6) in Table 



|are 






ii'i) 


= n(V'2) 


= Ai 




= n(V'4) 


= A2 




= n{^e) 


= A3 



a + 6 + c + V(a - 6 + c)2 + 8d2 
1 = 7, ~ l-U, 



a + 6 + c- v^(a - 6 + c)2 + Sd^ 



0.0, 



(7-6) 



7.4. Parameter set (4) 

In the parameter set (4), the G{XAg, 1, 1)3 state is most stable. The calculated 
values of o, b, c, d and / in (I6T7P are 

a = 0.75, 6 = 0.5, c = -0.1802, d = -0.1408, / = 0.3416. (7-7) 
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Then the occupation numbers for general spin orbitals V'j (j = Ij""" > 6) in Table 
IXTVl are 

n(V^i) = n(V'2) = Ai = a + Vc^ + d^ = 0.9786 ^ 1.0, 
n('03) = n(V'4) = A2 = a - v^c^Td^ = 0.521 « 0.5, 
n(^5) = A3 = 6 + \/2/ = 0.017 w 0.0, 

n(V'6) = A4 = 6 - V2f = 0.983 s:i 1.0. (7-8) 

The magnitude of spin density \s{m)\ is c + d + f = -0.1802 - 0.1408 + 0.3416 = 
0.0206. 

§8. Summary and discussion 

We applied the group theoretical bifurcation theory of the HF equation to the 
t2g-Hubbard model on a two-dimensional square lattice. By enumerating the axial 
isotropy subgroups of the R-reps of the group Go in the HF Hamiltonian space Whf 
in the cases of ordering vectors Qi,i = 0,1,2,3, we have obtained many types of 
broken symmetry states which bifurcate from the normal state through one step 
transition. 

It is shown that these states have various types of local order parameters (LOP): 
spin density s = (s^, s^, s^), quadrupole moment Q = {Q2, Q12, Q23, Q31), orbital an- 
gular momentum I = {li,l2,h), spin quadrupole moment = (Ql^, Qi2) Q23> Q31) 
and spin orbital angular momentum Z'^ = (^1,^25^3)' where A = 1,2,3. We have 
illustrated how the isotropy subgroup of a state can determine the canonical form 
of the HF Hamiltonian, the standard forms of the occupied spin orbitals and their 
occupation number, and types of LOP. 

We performed numerical calculations for all states discussed in §5 and §6 with 
some parameter sets. We found that various types of broken symmetric solutions 
can be the most stable state depending on parameter sets. Through these calcula- 
tions we found that non-collinear magnetic states:G{MEg, 1, 1)3, G^XB^g, 1, 1)3 and 
G{XAg, 1, 1)3 can be the ground state for some parameter sets. 

In §5 and §6 we considered broken symmetry states bifurcating through a single 
phase transition from the normal state. However as shown in §7 of I, there are many 
other states derived through two step transitions from the normal state. Here we 
consider two states bifurcatig from the G{MB2g, 0,0) treated in the Example 5.2. 
The isotropy subgroup of the G{MB2g, 0, 0) is 

G{MB2g, 0,0) = {E + T{ei)C2,)D2aLiSR (8-1) 

By a similar method to that in §7 of the previous paper, we obtain the following two 
types of isotropic subgroups 

G{MB2g,0, 0)a = {E + tu2y){E + T{ei)G2.)D2ahLiA{e3), 

G{MB2g, 0, 0)b = {E + tu2y){E + T{ei)C2xU2x)D2ahLiA{e^). (8-2) 
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We consider these states seperately. 

(a) G(M52<;,0,0)a state. 

This state has following non zero R}'^ 



00 



03 




R 



10 



R 



13 




(8-3) 



Thus from p-43p we obtain for m such that T{m) G Li 

Qi2{m) = 4I2C, Quim + Sj) = -Ahc, 
s^{m) = 2d + f, s^{m + ej) = 2d + f, 
Ql2{m) = 2126, Ql2{m + ej) = -2/26, 



(8-4) 



where j = 1,2. Thus we see that through the ferromagnetic transition of the 

G{MB2g,0, 0) state which has the anti-ferro quadrupole moment Q12, there appears 

anti-ferro spin quadrupole moment ^ ^^11 as ferro spin density s^. 

(b) G{MB2g,0,0)b state. 

This state has following non zero R''^ 



r! 
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13 




(8-5) 



e \ 
e \ R 
0// 

Thus we obtain, for m such that T{m) € Li and j = 1, 2, 

Qi2{m) = 4/2C, Qi2{m + ej) = -4/2C, 

s^{m) = 2d + f, s^{m + ej) = -{2d + f) 

Ql2im) = 2/26, Ql2im + ej) = 2/26. 

Thus we see that through the anti- ferromagnetic transition of the G{MB2g, 0,0) 
state, there appears ferro spin quadrupole moment (5x2 s>s well as anti-ferro spin 
density s^. 

By similar manner to that of above cases, we can see that states, which are 
derived through two step transition from the normal state, are coexisting states of 
{spin sensity , quadrupole moment and spin quadrupole moment }, or {spin sensity, 
orbital angular momentum and spin orbital angular momentum}. 
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Appendix A 

Derivation of quadrupole moment 

Here we derive ()3-36p for the quadrupole moment at the site m. The quadrupole 
moment operator Qij{m) at a site m is defined by^^^ 



"mps^mqs' 



3 2 

Qij{m) = ^ {4'mps\{'iXiXj - 5ijr'^)\4)mqs')a\r. 

p,q=l s,s' = l 
2 3 

s=l p,<?=l 

where (t)mps = (t>p{r - m)xs, Xs are the spin function: xi = I T),X2 = I i), and Qij 
is a 3 X 3 symmetric matrix whose {p, q) component is defined by 

i^i3)pq = j dr<l)p{r)i3xix, - r^6i,)cl)gir). (A-2) 
Thus the expectation value of the quadrupole moment at a site m given by 

2 3 

(A-3) 

s=l p,q=l 

The explicit forms of Qij are given by 

Qii = I -h I , Qi2 = Q21 



2/1 
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2/1 



Q22 = I 2/1 I , Q23 = Q32 

Q33 = I -h I , Q31 = = I I , (A-4) 

where 
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h = j dr(pi{r){x'^ - y'^)(j)i{r), 

l2 = 3 j dr(j)i{r)xy(j)2{r). (A-5) 

Thus from (jA-3p and ()A-4p we obtain 
2 

Qii(m) = /i ^ (2(aJ„i,a^iy) - (a^2sam2s') - (aL3s«m3s')) 

s,s'=l 
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2 



s,s'=l 
2 

Qssim) =h [-{(^LlsO-mls') - (aL2sam,2s') + 2(aL3s«rri3s')) ^^ss', 

s,s'=l 

2 

Ql2{m) = Q2l{m) = h Y ((aLlsam2.') + («L2s«rrils')) 

s,s'=l 
2 



s,s'=l 
2 

s,s'=l 

The spin quadrupole moment is defined by 

3 2 



Qiji-m) = \j2J2 i^rrvpsli^XiXj - Sijr'^)a^,,\(t)mqs')al^psa^g,, . (A-7) 

p,q=l s,s'=l 

Then we obtain (j3-38p for the expectation values of spin quadrupole moment. 



Appendix B 

Derivation of the orbital angular momentum 

Here we derive p-40p for the orbital angular momentum at the site m. Let 
Ij (j = 1,2,3) be the orbital angular momentum. Then operators of the orbital 
angular mormentum at the site m are expressed by 

3 2 

k{m) = J2 Yl (^■'nps\ii\(t>mqs')alnpsamqs' 
p,q=l s,s' = l 
3 2 

~ Yl {(pmps I ij I (Pmqs) O'lnps^mqs 
p,q=l s=l 
3 2 

~ y ^ y ']{'^i)pq'^rn,ps'^mqs- 0^'^) 
p,q=l s=l 

where 

{Li)pq = i'Pmpslhl'Pmqs) ■ (B"2) 

We denote a 3 x 3 matrix, whose {p,q) component is {Li)pg, by Lj. Within the t2g 
subspace we have^^^ 

000\ / -i \ / i 

Li = I i , L2 = , L3 = -i I . (B-3) 
-i \ i VOOO 
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Thus form (IB-ip and ()B-3|1 the expectation values li{m) of the orbital angular mo- 
mentum li{m) are given by 

ss' 
ss' 

h{m) = ^ii){{alnis'^m2s') - {<4n2s'^mis'))(^ss' ■ (B-4) 



The spin orbital angular momentum is defined by 

3 2 

i>mps I h/^ss' 1 4'mqs' ) '^Inps'^mqs' ■ 



iUm) 



p,q=l s,s'=l 

Then we obtain ()3-4ip for the spin orbital angular momenta. 



(B-5) 
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